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ABSTRACT 


",-\A  neutrally  stable  atmospheric  surface  layer  was 
suitably  simulated  In  a  low  speed  wind  tunnel  by  tripping 
the  boundary  layer  with  a  fence  and  letting  the  turbulent 
flow  develop  over  a  length  of  roughness  elements.  Jet 
exhaust  dispersion  characteristics,  simulated  by  a 
burner/no :zle  system,  were  Investigated  by  measuring  the 
horizontal  and  vertical  temperature  profiles  at  axial 
stations  downwind  from  the  nozzle  exit.  Dispersion 
sensitivity  to  different  nozzle  exit  conditions,  angles  cf 
Incidence  to  the  wind,  and  nozzle  surface  blockage  were 
Investigated.  The  results  were  compared  to  dispersion 
methods  used  in  the  Air  Quality  Assessment  Model  (AQAM). 
was  found  that  the  experimental  Jet  penetration  length  was 
much  shorter  than  that  assumed  in  AQAM,  and  that  a  plume 
rise  existed,  which  Is  not  Included  in  the  current  A  (JAM 
model.  Required  Inputs  of  the  Initial  dispersion 
coefficients  were  determined  as  a  function  of  wind 
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I.  INTRODUCTION 


The  Air  Force  has  developed  an  Air  Quality  Assessment 
Model  (AQAM)  to  predict  the  dispersion  of  pollutants  emitted 
from  typical  Air  Base  related  sources  (Ref.  1,  2  and  3). 

That  model  has  also  been  modified  to  better  simulate 
operations  from  a  Naval  Air  Station  (Kef.  4  ard  5).  The 
full  model  consists  cf  three  computer  programs:  Source 
Inventory,  Short  Term  and  Long  Term  Dispersion.  In  the 
Short  Term  Model  hourly  averaged  concentrat  ions  are 
predicted  over  a  grid  0f  receptors.  Air-base,  aircraft,  and 
off  air-base  (environ)  sources  are  considered  and  are 
modeled  as  point,  area,  or  line  sources  as  appropriate. 

In  the  Short  Term  Model,  the  Jet  exhausts  from  aircraft 
during  taxi  and  taiteoff  are  treated  as  line  sources  with  no 
plume  rise.  Initial  plume  dimensions  must  be  estimated  as 
well  as  the  distance  the  Jet  ’penetrates"  the  atmosphere 
before  coming  essentially  to  rest  relative  to  the  ambient 
wind.  Stationary  sources  such  as  engine  run  up  test  stands 
and  test  cells  are  treated  as  points  sources.  Likewise,  the 
initial  plume  dimensions  must  be  estimated  and  a 
determination  whether  plume  rise  is  a  factor  must  be  made. 
The  model  then  uses  these  input  data  to  determine  an  initial 
dispersion  coefficient  from  which  a  pseudo  upwind  point  or 


14 


line  source  Is  located.  The  sigma's,  or  dispersion 
coefficients,  are  based  on  the  stability  class  of  the 
atmosphere  and  the  downwind  distance  or  travel  time  for  an 
elevated,  relatively  low  velocity  source.  There  is  a 
question  then,  as  to  the  applicability  of  these  dispersion 
coefficients  to  turbojet  exhausts  which  are  often  sonic, 
have  high  shear  stress/turbulent  mixing  characteristics  and 
have  different  turbulence  scales  than  typical  exhaust 
stacks.  Evidence  that  there  is  a  discrepancy  has  been 
indicated  by  an  Air  1’orce  study  (Ref.  6)  which  presented  a 
pictorial  investigation  of  the  effects  of  atmospheric 
stability  on  a  Jet  exhaust.  It  is  evident  that  qualitative 
and  quantitative  data  are  needed  to  characterize  the 
dispersion  of  a  turbojet  exhaust  during  ground  operations. 
This  investigation  attempted  to  determine  turbojet  exhaust 
dispersion  as  a  function  of  Jet  characteristics  and  wind 
direction  under  laboratory  simulated,  neutrally  stable 
atmospheric  conditions. 


The  first  part  of  the  study  Involved  the  cor 
simulation  of  the  lower  portion  of  the  atmospher 
aeronautical  type  wind  tunnel.  Published  studie 
8  and  S)  have  indicated  that  it  is  not  sufficien 
simulate  only  the  mean  velocity  profile  of  the  a 
The  turbulent  structure  of  the  atmosphere  must  a 
correctly  modeled  simultaneously  with  the  mean  v 


rect 
e  in  an 
s  (Refs.  7, 
t  to 

tmosphere. 
Iso  be 
elocity 
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profile.  Proper  simulation  can  be  accomplished  following 
the  work  of  Sundaram  and  Ludwig  (Ref.  7),  tripping  the 
boundary  layer  with  a  suitable  fence  at  the  entrance  of  the 
wind  tunnel  test  section  and  then  letting  the  turbulent  flow 
develop  over  a  length  of  roughness  elements.  Suitable 
matching  between  the  boundary  layer  trip  and  the  roughness 
elements  then  results  in  an  area  of  constant  shear  stress 
which  adequately  models  the  atmospheric  flow. 

Once  an  atmospheric  surface  layer  has  been  suitably 
simulated,  a  turbojet  exhaust  must  then  be  simulated.  If  it 
is  assumed  that  heat  and  mass  have  the  same  turbulent 
diffusion  characteristics,  then  the  dispersion  rates  of 
pollutants  can  be  determined  by  measuring  the  temperature 
distributions  in  the  vertical  and  horizontal  directions. 

In  this  investigation  data  were  obtained  for  varying 
angles  of  incidence  to  the  oncoming  wind  for  choked  and 
unchoked  nozzle  flow  conditions,  for  two  Jet  temperatures, 
and  for  two  conditions  of  nozzle  surface  blockage.  The  data 
were  compared  to  the  dispersion  routines  used  in  the  AQAM 
computer  program. 
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II.  ATMOSPHERIC  SIMULATION 


In  order  to  make  meaningful  wind  tunnel  measurements  to 
determine  Jet  aircraft  exhaust  dispersion  rates  during  taxi, 
idle  and  takeoff  modes  of  operation  it  was  necessary  to 
correctly  model  the  lower  structure  of  the  atmosphere. 
Sundaram  and  Ludwig  (Ref.  7)  have  published  a  theoretical 
study  of  the  requirements  for  the  modeled  flow  and  an 
experimental  study  of  the  flow  generated  by  wind  tunnel 
techniques.  This  reference  was  used  as  the  principal  guide 
in  modeling  the  atmospheric  surface  layer  in  a  wind  tunnel 
at  the  Naval  Postgraduate  school.  Supplementary  discussions 
are  presented  in  Refs.  10-12 

A.  THEORETICAL  REQUIREMENTS  FOR  ATMOSPHERIC  MODELING 

Reference  13  describes  the  atmospheric  boundary  layer  as 
a  turbulent  layer  which  is  influenced  by  a  combination  of 
actions  due  to  surface  friction,  thermal  stratification,  and 
Coriolis  force.  The  atmospheric  boundary  layer  is  divided 
into  two  distinct  layers.  The  lower  region  called  the 
surface  boundary  layer  extends  up  to  a  nominal  height  of  100 
meters.  The  region  above  the  surface  layer  is  called  the 
Ekman  layer  and  extends  up  to  a  height  of  about  1000  meters. 
In  the  Ekman  layer  the  wind  direction  changes  with  height  so 
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that  an  equilibrium  balance  Is  maintained  between  the 
Coriolis,  pressure,  and  shear  forces.  In  the  surface  layer, 
however,  the  wind  direction  is  generally  considered  to 
remain  constant,  the  Coriolis  force  Is  not  an  Important 
parameter  and,  therefore,  the  variations  in  shear  stress  and 
vertical  momentum  with  Increasing  height  are  considered  to 
be  negligible.  The  Monin-Obukhov  theory,  Ref. 14,  and  other 
surface  layer  theories  make  the  assumption  that  the  flow 
characteristics  close  to  the  ground  can  be  expressed 
entirely  in  terms  of  conditions  at  the  surface.  This 
assumption  means  that  the  conditions  close  to  the  surface 
are  independent  of  the  conditions  In  the  Skman  layer.  The 
boundary  layer  equation  for  a  flat  plate  with  no  pressure 
gradient  is 


uau/3x  ♦  w9u/  3z  *  3t/3z. 


(1) 


Near  the  surface,  the  left  hand  side  of  equation  (1)  will  be 
small.  The  shear  stress  can,  therefore,  be  considered  to 
be  independent  of  height  in  the  region  close  to  the  surface. 
This  assumption  is  important  because  a  simple  theoretical 
description  of  a  neutrally  stable  atmospheric  boundary  layer 
is  made  possible  which  permits  the  simulation  of  such  a  flow 
in  an  experimental  wind  tunnel.  A  turbulent,  neutrally 


18 


stable  flow  structure  In  the  atmospheric  surface  boundary 
layer  can  be  described  entirely  by  the  following  parameters; 
the  kinematic  viscosity  v  ,  the  friction  velocity  u,  ,  arl 
the  roughness  parameter  zQ  (Ref.  7).  u,  is  defined  as  being 
equal  to  (  t©/ p  and  Zois  a  length  scale  which  describes 

the  Influence  of  a  rough  surface  on  a  particular  flow.  Most 
atmospheric  boundary  layers  are  observed  to  have  no  laminar 
sublayer  near  the  surface  and  are  described  to  be  fully 
aerodynamlcally  rough.  Using  the  above  parameters,  an 
aerodynamlcally  rough  flow  exists  if  (Ref.  7) 


Zq u,  /  v  >  3 .  (2) 

If  z0  is  <<  than  the  vertical  dimension  of  the  boundary 
layer,  then  the  mean  velocity  profile  of  the  wind  Is  given 
in  logrlthmlc  form  by  the  relationship  (Ref.  7) 


u/ur  «  ( l/<  )  In  (  z/Zq  )  ,  (3) 

where  u  is  the  average  velocity,  z  Is  the  vertical  distance 
from  the  surface  and  <  Is  Von  Karman's  constant.  This 
equation  Is  a  form  of  the  well  known  law  of  the  wall  and  is 
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Indicative  of  the  interaction  between  the  air  flow  above  a 
given  surface  and  that  surface.  The  existence  of  such  a 
velocity  profile  implies  (Ref.  7)  that  the  flow  is 
horizontally  homogeneous.  By  neglecting  the  Coriolis  force, 
a  neutrally  stable,  horizontally  homogeneous  flow  obeys  the 
following  relationship  (Ref.  7) 


3x/3z  *  (-1/  o  )  3  p/3x 


(4) 


which  describes  the  vertical  variation  of  shear  stress  with 
the  horizontal  pressure  gradient.  If  constant  shear  stress 
is  a  requirement  for  a  neutrally  stable  atmosphere,  then 
equations  (2)  through  (4)  identify  additional  requirements 
for  a  such  a  flow.  The  additional  requirements  are;  a  fully 
aerodynamlcally  rough  flow,  horizontal  homogeneity,  and  the 
absence  of  a  pressure  gradient.  To  accurately  model  the 
atmosphere  these  requirements  must  be  strictly  adhered  tc. 

If  these  requirements  are  strictly  maintained  In  a 
laboratory  flow,  then  it  can  be  expected  that  the  flow  will 
be  similar  to  the  atmospheric  surface  layer.  Strict 
horizontal  homogeneity,  however,  is  extremely  difficult  to 
obtain  in  the  laboratory.  Partial  simulation,  therefore, 
generally  has  to  be  imposed.  A  complete  discrlption  of  the 
effects  of  partial  simulation  is  given  in  reference  7  and 


will  not  be  repeated  here.  The  conclusion  reached  from 
reference  7  regarding  partial  simulation  was  that  a 
laboratory  flow  which  was  fully  aerodynamical ly  rough, 
nearly  horizontally  homogeneous  and  relatively  free  from 
pressure  gradients  was  a  reasonable  model  for  a  neutrally 
stable  atmosphere.  The  generation  of  such  a  flow  is 
described  in  the  next  section  of  this  report. 


3.  EXPERIMENTAL  REQUIREMENTS  TO  SIMULATE  ATMOSPHERIC  FLOWS 

In  the  previous  section  the  requirements  for  simulating 
neutrally  stable  atmospheric  flows  were  discussed.  That  is 
to  correctly  model  the  flow,  it  must  be  fully 
aerodynamically  rough,  horizontally  homogeneous,  free  of 
pressure  gradients  and  the  shear  stress  must  be  constant  in 
the  vertical  direction.  Aircraft  Jet  exhaust  dispersion  in 
the  "near  ground”  environment  was  the  primary  concern  of 
this  investigation*  therefore  only  the  surface  layer  cf  the 
atmosphere  was  modeled.  As  discussed  previously,  this  was 
possible  since  both  the  atmospheric  flow  and  the 
experimental  flow  are  described  completely  in  terms  cf 
conditions  at  the  surface  (ie.  the  parameters  u,  ,  Zq  , 
and  v  ).  The  similarity  requirements  were,  therefore,  not 
function  of  the  geostrophic  wind  in  the  real  flow  nor  a 
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function  of  the  free  stream  velocity  in  the  modeled  flow. 
Horizontal  homogeneity  cannot  he  strictly  satisfied  since 
for  an  experimental  flow  over  a  flat  plate,  equilibrium 
conditions  change  continually  as  a  function  of  downstream 
distance  (Ref.  7).  It  was  not  sufficient,  therefore,  to 
simply  compare  the  mean  velocity  profile  of  an  atmospheric 
flow  with  that  of  a  simulated  flow  at  one  location  within 
the  wind  tunnel.  The  additional  measurements  which  were 
required  for  partial  simulation  are  described  in  reference 
7.  Since  the  atmosphere  is  assumed  to  exhibit  a  logrithmic 
mean  velocity  profile,  this  part  of  the  similarity 
requirement  can  be  easily  verified.  In  reference  7  the 
turbulence  intensity  and  the  integral  scale  of  turbulence 
were  also  measured  and  compared  to  available  atmospheric 
data.  The  turbulence  intensity  measurements  in  this 
experiment  were  compared  to  the  results  of  reference  7. 
Integral  scales  of  turbulence  were  not  measured.  For 
virtually  Identical  experimental  set-ups  it  was  assumed  that 
if  all  other  experimental  measurements  compared  favorable  to 
those  of  reference  7,  then  the  scales  of  turbulence  would 
also  compare  favorably. 
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C.  EXPERIMENTAL  APPARATUS  AND  MEASUREMENTS 


1.  Wind  Tunnel  Modi  flea  ti ons 


The  experiment  was  conducted  in  a  lew  speed  wind 
tunnel  at  the  Naval  Postgraduate  School  (NFS),  Monterey,  Ca. 
The  wind  tunnel  draws  ambient  air  from  the  surrounding  area 
into  a  5  X  5  foot  test  section  approximately  21  feet  long. 

In  order  to  he  comparable  with  the  wind  tunnel  used  in 
reference  7,  the  wind  tunnel  at  the  Naval  Postgraduate 
School  was  modified  in  the  following  manner.  A  false 
ceiling  was  installed  to  maintain  an  approximate  zero 
pressure  gradient  in  the  axial  direction.  A  two  inch  high 
boundary  layer  trip  constructed  of  wood  was  placed  across 
the  entrance  to  the  test  section.  The  trip  was  followed  by 
2.5  X  2.5  feet  square  sections  of  roughness  elements  which 
covered  the  entire  floor  of  the  test  section.  The  roughness 

i 

elements  were  4  sided  regular  pyramids  which  were  3.75 
Inches  square  at  the  base  and  *).75  inches  in  height.  They 
were  constructed  out  of  1.5  inch  thick  fiberboard.  The 
tunnel  and  roughness  elements  are  shown  in  figure  1. 

2.  Pressure  Measurements 

Six,  1/S  Inch  diameter  pressure  taps  spaced  4  feet 
apart,  were  mounted  flush  along  the  centerline  of  the  wind 


tunnel  wall.  The  differential  area  required  to  naintain  a 
zero  pressure  gradient  in  the  axial  direction  was  determined 
from  an  initial  set  of  pressure  readings.  The  false  ceiling 
which  implemented  the  required  area  change  was  then 
installed  and  the  pressure  readings  were  again  taken  to 
verify  the  required  zero  axial  pressure  variation. 


3.  Electronic  Equicement 


All  velocity  related  measurements  were  taken  with  a 
Thermo  Systems  Inc.  linearized  hotwire  anemometer  system. 
DISA  single  wire  and  cross  wire  probes  and  probe  holders 
were  used  with  the  anemometer  system  (fig.  2  and  3).  A  DC 
digital  voltmeter  and  a  true  RMS  meter  were  used  to  record 
the  steady  and  fluctuating  single  wire  probe  data 
respectfully.  A  model  1015C  Thermo  Systems  Inc.  correlator 
was  used  in  conjunction  with  the  cross  wire  outputs  to 
obtain  the  sum  and  difference  of  the  two  signals  (fig.  4). 
The  true  RMS  meter  was  also  used  to  record  the  cross  wire 
outputs  via  the  correlator.  The  hotwire  probes  were  mounted 
under  the  wind  tunnel  such  that  the  probes  traversed  the 
wind  tunnel  in  a  vertical  direction.  A  Spectral  Dynamics, 
Spectra-Scope  Model  SD-330  frequency  counter  was  used  in 
conjunction  with  the  measurement  of  the  centerline  velocity 
(fig.  5). 
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4.  Wind  Tunnel  Centerline  Velocity  Measurements 


The  maximum  velocity  expected  in  the  wind  tunnel  was 
between  43  and  53  feet  per  second.  Existing  pitot-static 
systems  at  NFS  were  not  sensitive  to  tnis  small  velocity. 
Therefore,  the  relationship  between  the  Strcuhal  Number  and 
the  Reynolds  Number  was  used  to  determine  the  centerline 
velocity  (Ref.  15).  The  Strouhal  number  relates  the 
velocity  of  a  flow  tc  the  frequency  at  which  vortices  are 
shed  behind  a  circular  cylinder  which  is  immersed  in  that 
flow. 


Sn(Strouhal  Number)  =  (freq.  X  diam.  of  cyl.)  /  velocity. 


A.  2.2?5  inch  diameter  stainless  steel  wire  was  attached  to  a 
collar  which  was  fitted  over  a  single  hotwire  probe  such 
that  the  hotwire  was  in  only  one  wake  (shed  vortex  sheet)  of 
the  cylinder  (fig.  6).  The  probe,  collar  and  hotwire 
combination  was  placed  perpendicular  to  flow  at  the 
centerline  of  the  tunnel.  With  the  tunnel  on,  the  hotwire 
system  was  sensitive  to  the  velocity  of  the  free  stream  and 
the  output  was  observed  and  recorded  on  the  frequency 
counter.  The  observed  frequency  of  shed  vcrticies  and  the 
unique  relationship  between  the  Strouhal  Number  and  the 
Reynolds  Number  were  used  to  iteratively  determine  the 
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velocity  of  the  free  stream 


5.  Velocity  Profile  Measurements 

The  mean  velocity  profile  measurements  were  made 
using  the  single  hotwire  apparatus.  The  hotwire  system  was 
calibrated  in  the  linear  mode  such  that  the  centerline 
velocity  was  made  equal  to  one  volt.  The  single  wire  was 
also  calibrated  for  flow  direction  sensitivity.  This  was 
accomplished  by  rotating  the  probe  umtil  a  maximun  reading 
was  observed  on  the  DC  digital  voltmeter.  It  was  then 
assumed  that  the  probe  was  perpendicular  to  the  flow.  All 
subsequent  measurements  were  referenced  to  this  direction. 
The  linearized  output  was  averaged  through  a  ten  second  time 
constant  before  being  displayed  on  the  DC  digital  voltmeter. 
The  velocity  profile  was  obtained  by  recording  the  DC  value 
of  the  hotwire  output  as  it  traversed  the  wind  tunnel  in  the 
vertical  direction.  Vertical  measurements  were  referenced 
to  the  top  of  the  roughness  elements  and  were  measured  to 
within  an  accuracy  of  3.020  Inches  with  respect  to  that 
reference.  Velocity  profiles  were  obtained  at  two  axial 
locations.  The  locations  were  72  inches  and  109  inches 
downwind  from  the  boundary  layer  trip  at  the  entrance  to  the 
test  section.  The  two  axial  measurements  were  necessary  to 
ensure  the  required  degree  of  uniformity  in  the  Jet 
dispersion  region. 


. 
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6.  Turbulence/Shear  Stress  Measurements 


The  cross  wire  system  described  earlier  was  used  to 
measure  the  turblence  intensities  and  shear  stresses 
directly.  The  cross  wire  was  calibrated  such  that  the 
linearized  output  for  each  wire  was  of  equal  sensitivity. 
This  was  accomplished  by  first  orienting  one  wire  normal  to 
the  flow  and  then  adjusting  its  sensitivity  to  the  proper 
level.  The  probe  was  then  re-oriented  so  that  the  second 
wire  was  normal  to  the  flow  and  its  sensitivity  was  then 
matched  to  that  of  the  first  wire.  For  data  acquisition, 
the  probe  was  oriented  axially  such  that  the  longitudinal 
turbulence  intensity  was  proportional  to  the  RMS  value  of 
the  sum  of  the  cross  wire  outputs.  Likewise,  for  the  same 
orientation,  the  vertical  turbulence  intensity  was 
proportional  to  the  RMS  value  of  the  difference  of  the  cross 
wire  outputs.  The  turbulent  shear  stress  is  defined  as 
-  ou'w’ .  This  shear  stress  was  obtained  directly  since  u'w' 
is  proportional  to  the  square  of  the  sum  minus  the  square  of 
the  difference  of  the  cross  wire  outputs.  Like  the  single 
wire  data,  the  cross  wire  data  was  recorded  as  the  probe  was 
traversed  in  the  vertical  direction  at  the  same  two  axial 
locations.  For  redundancy,  the  vertical  profile  of 
turbulence  intensity  was  also  measured  using  the  single  wire 
probe. 
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D.  RESULTS  AND  DISCUSSION 


All  the  measurements  obtained  at  the  two  axial 
locations,  le.  72  and  1J9  Inches  from  the  entrance  of  the 
wind  tunnel  test  section,  were  taken  with  a  wind  tunnel 
center  line  velocity  of  approximately  4k}  feet  per  second. 

1.  u,  and  zQ 

From  the  law  of  the  wall, 

u/u,  =  5.75  log  z/zQ  (5) 


Equation  (5)  can  be  rewrlten  as 


u/Uco  «  (u,  /Uoo  )  5.75  ( log  z  -  log  z0  )  .  (6) 


When  u  /Uoo  Is  plotted  vs.  log  z, 


slope  »  5.75  (u,/Uoo  ). 


i  7 ) 
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u/Uco  TS.  log  z  was  plotted  for  both  axial  positions  and  a 
least  squares  fit  was  obtained  through  the  data  points 
(figures  7  and  S).  The  slopes  were  measured  and  u,  was 
determined  by; 


u,  =  (slope  •  U®  )/5.?5. 


(a) 


u,  values  of  2.32  and  2.52  feet  per  second  were  obtained  at 
the  72  and  109  inch  positions  repectively.  Using  the  least 
squares  fit  data,  the  calculated  u,  was  used  to  determine  zQ 
by; 


log  z0  *log  z  -  u/ { 5 . 75  •  u,  )  (9) 

z0  values  of  0.0215  and  0.0226  were  obtained  at  72  and  109 
inches  respectively. 

2*  Vertical  Distribution  of  Shear  Stress 


Since  Up  is  proportional  to  the  shear  stress. 


a  measure  of  the  average  shear  stress  for  each  radial 
position  can  be  obtained  from  the  mean  velocity  profile 
data.  Non-dimensionalizlng  by  \Jco  » 


u'w'  /  U®1  *  (u.  /U»  )* 


(11) 


The  measured  shear  stress,  u'w'»  from  the  cross  wire  data, 
was  also  non-dimentlonalized  by  \Jao  and  plotted  as  a 
function  of  the  vertical  distance  above  the  roughness 
elements  for  both  axial  positions.  Figure  9  presents  the 
cross-wire  data.  Also  shewn  are  the  values  of  average 
stress  determined  from  the  law  of  the  wall. 

Figures  10  and  11  compare  the  results  with  those 
presented  in  reference  7.  It  is  observed  that  the  shear 
stress  continued  to  increase  until  very  near  the  wall 
whereas  the  shear  stress  from  reference  7  began  to  decay 
further  from  the  wall.  Figure  10  shows  that  the  shear 
stress  was  approximately  constant  for  heights  from  2  to  5 
inches  above  the  roughness  elements  at  the  72  inch  position. 
At  the  109  inch  position  (fig.  11),  the  shear  stress  decayed 
gradually  from  the  near  wall  region.  In  an  attempt  to 


modify  the  near  wall  shear  stress  distribution,  different 
combinations  of  trip  heights  and  wind  speeds  were  tried. 

This  proved  to  be  unsuccessful.  The  mean  velocity  profile 
was  observed  to  vary  much  more  rapidly  in  the  axial 
direction  and  the  relatively  constant  shear  stress  region 
obtained  with  the  two  Inch  trip  was  destroyed  in  varying 
degrees  depending  on  which  trip  or  what  tunnel  speed  was 
used.  There  was  also  a  somewhat  greater  variation  in 
horizontal  homogeneity  In  this  experiment  relative  to 
reference  7.  This  was  reflected  in  the  increased  slope  of 
figure  S  over  that  of  figure  7.  This  was  also  reflected  In 
figure  9  by  noting  that  the  average  shear  stress  differed  by 
about  20  percent  between  the  two  positions.  Moreover,  :Q 
Increased  by  about  50  percent  over  the  axial  range  compared 
to  about  30  percent  for  the  experiment  of  reference  7.  It 
was  also  noted  In  reference  7,  however,  that  the  degree  of 
uncertainty  in  the  atmosphere  may  be  much  larger  than  the  30 
percent  figure. 

3.  Vertical  and  Horizontal  Turbulence  Intensities 

Figures  12  and  13  compare  the  measured  vertical  and 
horizontal  turbulence  intensity  measurements  at  the  two 
axial  locations  with  the  data  of  reference  7.  In  general, 
the  vertical  and  horizontal  data  agreed  extremely  well  with 
the  data  of  reference  7.  The  horizontal  turbulence  was  a 
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little  lower  for  both  positions  in  the  range  of  z/zQ  from  1 
to  30  or  40.  This  corressponded  to  a  height  of  about  one 
inch  above  the  roughness  elements  and  was  the  area  of 
difference  noted  in  figures  10  and  11. 

Since,  (1)  a  region  of  constant  shear  stress  was 
obtained  in  good  agreement  with  reference  7,  (2)  the 
turbulence  intensities  were  relatively  constant  over  the 
same  region,  (3)  the  roughness  parameter  z0  was  comparable 
to  that  of  reference  7,  and  (4)  equation  (2)  was  satisfied, 
the  simulated  atmospheric  surface  layer  was  considered  to  be 
adequate  for  the  initial  study  of  Jet  exhaust  dispersion. 
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Figure  2 

Single  Hotwire  probe 


Figure  3 
Cross  Hotwire 
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Figure  4 

Anemometer  and  Correlator 


Figure  5 

Frequency  Counter 
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Velocity  Profile  at  72  inches  Downwind 


Figure  8 

Velocity  Profile  at  109  inches  Downwind 
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Figure  9 

Vertical  Distribution  of  Shear  Stress  at 
72  and  109  inches  Downwind 
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Figure  10 

Vertical  Distribution  of  Shear  Streaa  at 
72  Inches  Downwind  Compared  to  the  Data  of  Ref.  7 
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Figure  11 

Vertical  Distribution  of  Shear  Stress  at 
109  Inches  Dovnvind  Compared  to  the  Data  of  Ref.  7 
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Figure  13 

Comparison  of  Turbulence  Intensities  at 
109  inches  to  the  Data  of  Ref.  7 


III.  SIMULATED  JET  DISPERSION  IN  A  WIND  TUNNEL 


A.  BACKGROUND 

The  AQAM  model  treats  most  aircraft  emissions  as  finite 
line  sources.  Aircraft  engine  test  cells  and  run  up  stands 
are  treated  as  point  sources.  Line  sources  are  modeled  by  a 
general  line  source  or  puff  type  model  (Ref.  1-3).  Point 
sources  are  modeled  by  the  application  of  Gaussian  plume 
theory.  This  theory  represents  the  downwind  concentration  of 
pollutants  from  an  elevated  source  as  Gaussian  distributions 
in  both  the  horizontal  and  vertical  directions.  When 
applied  to  aircraft  exhausts,  the  line  source  method  first 
assumes  that  the  emissions  from  the  moving  aircraft  come  to 
rest  relative  to  the  ambient  wind  after  a  penetration  length 
due  to  the  velocity  effects  of  the  Jet.  After  coming  to 
rest,  the  line  source  is  generally  -segmented .  The  segmented 
portions  of  the  line  source  are  then  assumed  to  behave  as 
pseudo-upwind  line  sources  which  are  dispersed  downwind  in  a 
Gaussian  manner.  Both  methods  of  modeling  use  dispersion 
coefficients  oy  and  ax  which  are  the  standard  deviation 
points  of  a  Gaussian  or  normal  distribution  curve  and  are 
described  by  Turner  in  his  workbook  on  dispersion 
estimations  (Ref.  16).  As  mentioned  earlier,  the  dispersion 
coefficients  are  based  on  the  atmospheric  stability  class 


and  the  time  or  distance  traveled  downwind  from  a  relatively 
low  velocity  source.  The  applicability  of  these 
coefficients  relative  to  a  high  velocity,  horizontally 
emitted  source,  such  as  a  turbojet  engine,  was  investigated 
for  various  Jet  parameters. 

B.  EQU IP5MENT 

1.  General  Discussion 

The  investigation  was  conducted  in  the  low  speed 
wind  tunnel  described  in  the  previous  section.  A 
burner/nozzle  system  provided  the  source  for  a  simulated  Jet 
exhaust.  The  temperature  profile  of  the  exhaust  was 
measured  since  it  was  assumed  that  mass  and  temperature  are 
dispersed  by  the  same  mechanism.  The  temperature  profiles 
were  obtained  in  the  vertical  and  horizontal  directions 
using  a  thermocouple  apparatus. 

2.  Burner/Nozzle 

Oxygen  and  ethylene  were  mixed  and  ignited  at  the 
entrance  to  a  burner  can,  figure  14.  The  mixture  was 
augmented  with  low  pressure  air  for  burner  cooling  and  to 
create  the  desired  flow  rates.  The  hot  gases  were 
discharged  through  a  1/8  inch  diameter  stainless  steel  tube 
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into  the  flow  within  the  wind  tunnel.  The  tube  was  inserted 
through  the  tunnel  floor.  It  was  made  with  one  90  degree 
bend  which  aligned  the  exhaust  parallel  to  the  tunnel  floor. 
A  nozzle  was  created  by  reducing  the  exit  diameter  of  the 
tubing  to  0.058  inches.  The  tube  was  positioned  3  inches 
above  the  roughness  elements.  This  position  was  chosen 
since  it  was  about  the  center  of  the  area  of  constant  shear 
stress  discussed  in  the  atmospheric  simulation  section  of 
this  report.  Stagnation  pressure  and  temperature  probes 
were  used  to  determine  the  nozzle  exit  conditions.  The 
nozzle  entered  the  wind  tunnel  at  a  position  72  inches 
downstream  from  the  boundary  layer  trip  described  in  the 
previous  section. 

3.  Temperature  Measuring 

Horizontal  and  vertical  temperature  profiles  were 
obtained  using  a  chromel-alumel  thermocouple  probe.  The 
thermocouple  was  referenced  to  an  Omega  electronic  Miniature 
Cold  Junction  (MCJ).  The  nozzle  and  thermocouple  system  is 
shown  in  figure  15.  A  probe  holder  and  a  traverse  system 
permitted  positioning  of  the  thermocouple  probe  in  varying 
axial,  transverse,  and  vertical  locations.  The  traverse 
system,  an  electrical-mechanical  device  mounted  under  the 
wind  tunnel  (fig.  16),  facilitated  positioning  in  the 
vertical  and  transverse  direction.  The  probe  holder, 
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mounted  on  the  traverse  system,  permitted  the  thermocouple 
probe  to  be  positioned  in  the  axial  direction.  Two  X-Y 
plotters  (fig.  17)  were  used  to  record  the  temperature 
profiles  as  a  function  of  downwind  distance.  One  plotter 
was  used  for  the  vertical  temperature  profile  and  the  other 
one  was  used  for  the  horizontal  temperature  profile.  Each 
plotter  was  calibrated  for  displacement  in  inches  and 
temperature  in  milll-volts. 

C.  DISCUSSION  OF  RESULTS 

1 .  Introduction 

The  dispersion  data  from  the  experimental  Jet 
eihaust  were  used  to  determine  1)  If  the  plume  dispersion 
was  Gaussian,  2)  the  sensitivity  of  the  Jet  dispersion  to 
the  wind  direction,  3)  the  sensitivity  of  the  Jet  exhaust 
dispersion  to  nozzle  exit  Mach  number  (choked  and  unchoked) 
under  no-wind  and  wind  conditions,  4)  the  sensitivity  of 
the  Jet  dispersion  to  upwind  nozzle  surface  blockage  and  5) 
the  effect  of  increased  jet  exhaust  exit  temperature  on 
plume  rise.  The  first  four  items  were  investigated  with  an 
exhaust  stagnation  temperature  at  550  degrees  F  while  the 
stagnation  temperature  for  the  plume  rise  investigation  was 
increased  to  700  degrees  F.  For  choked  flow  the  exhaust 
stagnation  pressure  was  32  psia.  For  unchoked  flow  it  was 
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22  psla. 


2.  Attainment  of  Gaussian  Dispersion 

A  Gaussian  or  normal  distribution  with  a  mean  value 
of  zero  is  defined  by  the  equation 


y  =  (1/  (JZT*0')}  •  exp[-.5(x/<rf  ]  . 


(12) 


A  standard  normal  distribution  (fig.  16)  is  a  normal 
distribution  in  which  <r  is  set  equal  to  one.  In  this  case, 
the  maximun  value  of  the  ordinate  of  equation  12  is  1  /  J2  rr 
or  0.3969.  The  value  of  the  ordinate  at  one  standard 
deviation,  where  i  =  cr  ,  of  equation  12  is  0.242.  To 
determine  whether  the  measured  temperature  profiles  were 
Gaussian  the  profiles  were  standardized  in  the  following 
manner.  For  any  radial  temperature  profile  the  peak 
temperature  (Tmax)  was  assumed  to  be  at  the  mean  of 
the  distribution.  That  temperature  was  scaled  such  that 
Tmax/C  was  equal  to  0.3989.  The  ordinate  values  of  the 
temperature  profile  were  then  scaled  by  the  factor  C.  The 
ordinate  value  of  the  one  point,  T  ^  ,  of  the  profile  was 
determined  by?  T <r  *  00.242.  At  this  temperature  the 
width  of  the  temperature  plume  is  equal  to  2  o’.  The 
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standard  deviation  which  was  determined  in  this  manner  was 
used  to  non-dimensionall ze  or  standardize  the  abscissa  of 
the  temperature  profile.  The  standard  deviations  are 
tabulated  in  Tables  I,  II  and  III  for  each  of  the  conditions 
of  the  experiment.  The  standardized  data  points  were  then 
plotted  against  a  standard  normal  distribution  curve  and 
compared  for  closeness  of  fit.  Figures  19  through  24 
compare  the  vertical  and  horizontal  temperature  profiles  for 
a  choked  nozzle  flow  to  the  standard  normal  distribution 
curve.  The  Jet  exhaust  was  oriented  at  zero  degrees 
incidence  to  the  free  stream  wind  direction.  The  wind 
tunnel  centerline  velocity  was  40  feet  per  second.  The 
different  sets  of  data  represent  measurements  taken  at 
different  axial  positions  downwind  from  the  Jet  exhaust. 

The  downwind  distance  was  expressed  in  terms  of  Jet 
diameters.  The  plume  was  near  Gaussian  after  approximately 
eight  Jet  diameters  downstream  from  the  exit  plane  of  the 
Jet  for  both  the  horizontal  and  vertical  profiles.  The 
distance  required  for  the  profile  to  become  Gaussian 
appeared  to  be  the  result  of  the  expansion  process  of  the 
fluid  as  it  exited  the  nozzle.  Similar  plots  were  obtained 
for  the  other  zero  incidence  conditions  shown  in  Tables  I 
and  II.  For  each  case  the  plume  exhibited  near  Gaussian 
properties  at  about  eight  Jet  diameters  downstream.  The 
dashed  line  on  Tables  I  and  II  indicate  the  position  at 
which  the  plume  exhibited  Gaussian  characteristics. 
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As  the  nozzle  was  rotated  to  the  wind  direction,  the 
Jet-wind  interactions  caused  the  horizontal  profile  to  be 
negatively  skewed.  The  degree  of  skewness  diminished  as  a 
function  of  downwind  distance  as  the  effect  of  the  jet  exit 
velocity  was  overcome  by  the  force  of  the  wind.  At  the 
point  where  the  wind  was  the  dominant  dispersion  factor  (ie. 
past  the  penetration  length)  the  plume  again  exhibited 
Gaussian  charac teri s tics .  ?lgure  25  compares  the  data 
obtained  at  an  axial  distance  of  1.72  jet  diameters  for 
different  angles  of  incidence  to  the  wind.  The  angles 
investigated  were  25,  45,  67  and  90  degrees  and  the  data 
were  ncn-dimensi onalized  in  the  same  manner  as  discussed 
previously.  The  degree  of  skewness  in  figure  25  is  evident 
fcr  each  angle  of  incidence  when  data  are  compared  to  the 
standard  normal  distribution  curve.  Similar  data  for 
different  downwind  stations  are  plotted  on  figures  26 
through  32.  The  vertical  temperature  profiles  exhibited  no 
skewness  and  were  very  nearly  Gaussian  at  each  data 
position.  The  o'  ' s  which  were  obtained  for  the  vertical  and 
horizontal  dispersions  are  tabulated  in  Table  III.  The 
downwind  distance  at  which  the  plume  became  Gaussian  is 
indicated  for  each  angle  by  a  dotted  line. 

3.  Sensitivity  of  Exhaust  Dispersion  to  Angle  of 
Incidence  to  the  wind. 
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A  spatial  representation  of  the  Jet  exhaust  as  a 
function  of  downwind  and  horizontal  distances  from  the 
centerline  of  the  nozzle  are  plotted  on  figure  33  for  jet 
exhaust  angles  of  incidence  to  the  wind  of  29,  45,  67  and  90 
degrees.  As  expected,  as  the  angle  of  incidence  was 
increased,  the  length  normal  to  the  wind  direction  where  the 
dispersion  distribution  was  essentially  uniform  (ie.  sigma 
approaches  infinity)  increased  (in  a  nonlinear  relationship, 
fig.  34).  The  downwind  distance  at  which  the  plume  became 
approximately  uniformly  distributed  was  extremely  sensitive 
to  the  angle  of  incidence.  The  downwind  distance  decreased 
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quite  rapidly  as  the  angle  of  incidence  was  increased  in  an 
almost  linear  relationship  (fig.  34).  From  Table  III  it  can 
also  be  seen  that  the  corresponding  sigmas  increased 
considerably.  Quantitative  data  was  attempted  for  angles  of 
incidence  greater  than  90  degrees.  The  probe  could  not  be 
translated  further  upstream  than  the  90  degree  Jet 
orientation.  At  this  position,  the  exhaust  of  the  jet  was 
completely  dispersed.  Qualitatively  then,  as  the  angles  of 
incidence  were  increased,  the  exit  jet  momentum  and  wind 
action  acted  together  to  rapidly  disperse  the  exhaust  in  the 
downwind  direction.  The  dispersion  coefficients  were 
increased  which  resulted  in  a  decreased  downwind  distance 
where  the  plume  became  uniformly  distributed. 


4 .  Sensitivity  of  Dispersion  Rate  to  Nozzle  Eiit  Mach 
dumber  (choked  and  unchoked),  for  Wind  ana  No-wind 
Conditions 


Table  I  presents  the  results  in  terms  of  the  size  of 
the  one  standard  deviation  point  of  the  plume  in  jet 
diameters  as  a  function  of  downwind  distance  from  the  nozzle 
exit.  It  can  be  seen  that  the  plume  dimensions  were 
relatively  insensitive  to  varying  exit  conditions  (choked, 
unchoked  and  high  temperature).  The  plume  dimensions  for 
both  the  horizontal  and  vertical  directions  varied  less  than 
one  Jet  diameter  between  each  condition  over  a  downwind 
distance  of  about  120  Jet  diameteres.  The  data  for  the 
unchoked  nozzle  with  wind  dispersed  about  the  same  as  did 
the  choked  flow.  For  choked  flow  the  Mach  number  at  the 
exit  was  1.0.  For  the  unchoked  flow  the  Mach  number  was 
about  0.8.  Therefore,  both  had  high  velocities  relative  to 
a  typical  stack  exhaust. 

The  dispersions  from  a  choked  nozzle  with  and 
without  wind  were  also  measured.  The  two  sets  of  data  for 
both  the  horizontal  and  vertical  directions  are  plotted  in 
figures  35  and  36.  It  is  seen  that  the  dispersion  of  the 
plume  became  affected  by  the  wind  structure  at  about  30  jet 
diameteres  downwind  from  the  nozzle  exit.  Between  the 
nozzle  exit  and  about  30  jet  diameters  downstream,  the 
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no-wind  and  wind  data  compare  fairly  well  with  each  other. 
Further  downwind  from  that  position,  however,  the  two  curves 
diverge.  It  was  apparent  then  that  the  exit  conditions  of 
the  nozzle  (ie.  the  penetration  length)  governed  the 
dispersion  of  the  plume  for  about  the  first  30  Jet  diameters 
downwind,  and  thereafter  the  dispersion  was  primarily  a 
function  of  the  environmental  conditions. 


5.  Effect  of  Nozzle  Surface  Blockage 


To  investigate  the  effect  of  flow  disturbances 
(caused  by  the  tubing  upwind  of  the  nozzle)  on  the 
dispersion  of  the  Jet  exhaust,  two  different  non-heat 
conducting  disks  were  separately  fitted  on  the  nozzle 
tubing.  These  blockage  devices  increased  the  upstream  area 
relative  to  the  1/6  inch  tubing  by  a  factor  of  15  and  75 
respectfully.  Table  II  contains  the  results  and  it  was 
observed  that  upstream  blockage  had  no  significant  effect  on 
the  dispersion  rates.  These  upstream  blockage  effects  were 
apparently  dissipated  by  the  time  that  the  Jet  had 
penetrated  approximately  30  Jet  diameters  aft  of  the  nozzle. 
These  data  Indicate  that  aircraft  configuration  may  not 
significantly  affect  the  Jet  dispersion  for  the  slowly 
moving  taxi  operation. 


As  stated  previously,  the  effect  of  exit  temperature 
had  little  effect  on  the  plume  dimensions  as  a  function  of 
downwind  distance.  An  investigation  of  the  plume  rise  was 
conducted  by  measuring  the  distance  between  the  peak 
temperature  point  of  the  temperature  profile  and  the 
centerline  of  the  nozzle  for  two  exhaust  temperatures.  The 
two  sets  of  data  are  compared  as  a  fuction  of  downwind 
distance  in  figure  37.  Because  of  the  small  distances  it 
was  difficult  to  make  accurate  quantitative  conclusions  with 
respect  to  the  plume  rise.  From  the  data,  however,  it 
appeared  that  the  plume  rise  was  about  one  Jet  diameter  per 
1  <90  degrees  F  change  in  exit  temperature.  Equally  as 
interesting  as  the  plume  rise  was  what  appeared  to  be  a 
small  degree  of  looping  caused  by  the  turbulent  structure  of 
the  atmosphere.  This  phenomenon  is  shown  in  figure  37  as  a 
dotted  line.  The  looping  phenomenon  is  also  described  in 
reference  6. 


D.  COMPARISON  OF  EXPERIMENTAL  RESULTS  WITH  AQAM 

The  line  source  dispersion  employed  for  aircraft  exhaust 
Jets  in  AQAM  first  assumes  that  the  Jet  exhaust  has  a 
penetration  length  equal  to  l-4«3  meters  (ie.  approximately 
140  Jet  diameters).  The  penetration  length  is  that  length 
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where  the  velocity  effects  of  the  exhaust  cease  to  be  a 
factor  in  dispersing  the  exhaust.  Thus,  the  exhaust 
essentially  comes  to  rest  with  respect  to  the  ambient  air 
mass.  At  this  distance  then,  the  exhaust  is  being  dispersed 
entirely  by  the  action  of  the  wind.  At  the  penetration 
length,  estimates  of  the  plume's  lateral  and  vertical 
dimensions  a*e  made.  The  lateral  width  w,  and  vertical 
height  A  z  then  define  the  dimensions  of  a  line  source  of 
lenth  L.  The  lateral  width  is  assumed  to  remain  constant 
over  the  length  of  the  line.  Once  the  line  source  has  been 
defined,  the  line  is  generally  segmented  and  dispersed 
downwind  using  the  Gaussian  theory  of  dispersion.  That  is» 
from  the  dimensions  cf  the  line  source  an  initial  and  alo 
are  determined  from  empirical  relationships  oyo=  w/2.4  and  ctzo 
=  A  z/2.4  respectively.  With  these  values  of  oy*  and  azo  an 
upwind  distance  for  a  pseudo  upwind  line  source  is 
determined  using  the  oy  and  at  vs.  distance  or  time  charts 
of  reference  7.  The  pollution  concentrations  are  then 
dispersed  from  the  determined  line  source  in  a  Gaussian 
manner.  If  w  and ^2  are  not  specified  inputs,  the  AQAM 
model  uses  default  values  of  22  meters  and  8  meters 
respectfully. 

The  penetration  length  determined  in  this  investigation 
for  zero  degrees  orientation  to  the  wind  was  approximately 
30  Jet  diameters,  considerably  shorter  than  the  140  Jet 
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diameters  of  AQAM.  The  30  Jet  diameter  figure  was  further 
investigated  using  the  dispersion  data  from  tests  with  the 
nozzle  oriented  at  other  than  zero  degrees  to  the  wind  (fig. 
38).  The  30  jet  diameter  penetration  distance  determined 
from  the  axial  conditions,  indicated  by  an  arc  in  figure  38, 
was  compared  to  the  point  where  the  dispersion  curves 
exhibited  no  skewness  (ie.  where  the  distributions  were 
Gaussian  normal  to  the  wind  direction).  It  was  assumed  that 
beyond  this  position  the  jet  effects  could  be  neglected  and 
that  the  Jet  would  be  dispersed  by  the  ambient  wind  only. 

This  assumption  is  not  entirely  correct  since  the  plume 
continued  to  exhibit  horizontal  displacements  normal  to  the 
wind.  Along  the  Jet  axis  the  penetration  length  increased 
as  a  function  of  incidence  to  the  wind  from  about  28  to  58 
Jet  diameters,  still  considerably  less  than  the  140  jet 
diameters  assumed  by  AQAM.  This  relationship  is  plotted  in 
figure  3S  and  appears  to  be  approximately  linear  between  29 
and  90  degrees  relative  to  the  axial  orientation.  The 
distance  downwind  from  the  nozzle  where  the  plume  became 
Gaussian,  however,  remained  approximately  constant  at  25-30  jet 
diameteres.  It  was  also  noted  that  the  corresponding  cr's  at 
the  Gaussian  point  for  each  angle  became  wider  with 
increased  angles  of  incidence  due  to  more  turbulent  mixing. 

This  relationship  is  plotted  in  figure  40.  With  the  limited 
number  of  data  points  available  it  appears  that  the 
relationship  follows  a  smooth  curve.  <ry  and  ai  and  the  Jet 
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penetration  length  seem  to  be  predictable  between  zero  and 
9d  degrees  of  incidence  to  the  wind  for  neutral  stability 
and  these  nozzle  exit  conditions.  These  values  could  be 
used  in  AQAM  to  determine  w  and  Az  instead  of  estimating 
them  or  using  the  default  values.  The  relationship  shown  in 
figure  39  could  also  be  used  as  input  for  Jet  penetration 
length.  However,  since  the  penetration  lengths  are  small, 
it  may  be  better  to  neglect  this  effect  entirely. 


A  comparison  of  the  rates  of  dispersion  for  zero  degrees 
Incidence  to  the  wind  (i.e.  the  slope  of  the  and  civs. 
distance  curves)  was  made  between  a  neutrally  stable 
atmosphere  (class  D)  in  reference  1  and  the  experimental 
data.  As  can  be  seen  from  figures  35  and  36  the 
experimental  dispersion  rates  were  greater  than  the  class  D 
rate  used  by  reference  1  for  a  neutrally  stable  atmosphere. 
Reference  1  indicated  that  this  may  be  expected.  The 
dispersion  rates  were  obtained  from  data  compiled  from 
typical  low  velocity  stacK  emissions  and  not  from  high 
velocity  sources.  In  general  terms,  it  appears  that  for  the 
conditions  of  this  experiment,  a  more  realistic 
representation  of  the  actual  dispersion  rate  could  be 
obtained  if  the  atmospheric  stability  class  was  decreased  by 
one . 


For  a  vertical  stack. 


AQAl“i  also  predicts  the  plume  rise 
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due  to  thermal  effects  by  one  of  two  relationships.  These 
relationships  were  Investigated  to  see  if  they  could  be  used 
to  predict  the  plume  rise  observed  in  the  experiment.  For 
the  HOLLAND  plume  rise  formula,  a  rise  of  3.05  Jet  diameters 
per  100  degrees  F  was  predicted  and  for  the  CARSON-fiOSES 
formula  a  rise  of  0.008  Jet  diameters  was  predicted  per  100 
degrees  F.  This  compares  to  about  1  Jet  diameter  plume  rise 
observed  in  the  experiment.  Again,  the  discrepancy  seems  to 
be  in  the  use  of  equations  developed  for  low  velocity 
sources  for  sources  with  high  velocity. 


Figure  15 

Exhaust  Nozzle  and  Thermocouple  Probe 


Figure  16 

Thermocouple  Traverse  System 
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Standard  Normal  Distribution 


Downwind  Distance 


Temperature  to  a  Standard  Normal  Distribution 
for  Zero  Degrees  Angle  of  Incidence 
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Temperature  to  a  Standard  Normal  Distribution 
for  Zero  Angle  of  Incidence 
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25.8  Jet  Diameters  Downwind  from  the 


of  Incidence 


86.2  Jet  Diameters  Downwind  from  the 


Uniform 


Spatial  Representation  of  Plume  Width  as  a  Function  of 
Downwind  and  Horizontal  Distances  from  the  end  of  the 
Nozzle  for  Different  Angles  of  Incidence  to  the  Wind 


Downwind  Distance 
Where  Plume  is 
Uniformly  Distributed 


No  wind  condition 
Turner  Stability  Class 


Rise  vs.  Downwind  Distance 


Penetration  vs.  Downwind  and  Horizontal  Distance  from 
Nozzle  for  Different  Angles  of  Incidence  to  the  Wind 


Angle  of  Incidence  (degrees) 


Figure  39 
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Downwind  Distance  JD 
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SLgma  Vert, 
(choked ) 
Blockage 

(JD) 

Sigma  Horz. 
(choked) 
Blockage 

(JD) 

none 

small 

large 

none 

small 

large 

1.72 

.574 

.574 

.574 

.862 

.862 

.862 

8.62 

1.00 

1.00 

1.00 

.862 

.862 

1.29 

17.3 

1.43 

1.43 

1.72 

1.72 

1.72 

1.72 

25.8 

2.30 

2.30 

2.30 

2.58 

2.60 

3.01 

34.5 

2.90 

3.16 

3.44 

3.01 

3.00 

3.87 

51.7 

4.02 

4.59 

4.59 

4.74 

4.90 

5.17 

86.2 

6.61 

7.04 

7.75 

7.32 

7.50 

7.75 

120.6 

9.19 

.19 

9.19 

11.20 

11.30 

10.34 

Table  II 


1-Sigma  Plume  Width  of  the  Vertical  and  Horizontal 
Temperature  Profiles  at  Zero  Degrees 
Incidence  to  the  Wind  for  Different 
Nozzle  Blockages 
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Sigma  Vert.  (JD)  II  Sigma  Horz.  (JD) 


1-Sigma  Plume  Width  of  the  Vertical  and  Horizontal 
Temperature  Profiles  for  Different  Angles  of 
Incidence  to  the  Wind 


IV.  SUMMARY  AND  RECOMMENDATIONS 

In  summary,  a  neutrally  stable  atmospheric  surface  layer 
was  suitably  simulated  in  a  low  speed  wind  tunnel. 

Horizontal  and  vertical  temperature  profiles  aft  of  a  nozzle 
which  was  immersed  in  that  flow  were  measured.  The 
dispersion  data  from  the  simulated  jet  exhaust  were 
evaluated  for  different  nozzle  conditions  and  compared  to 
the  dispersion  data  used  in  the  AQAM  Model. 

•  It  is  recommended  that  additional  studies  be  made  for 

conditions  of  other  than  neutral  stability.  Measurements 
made  in  less  stable  surface  layers  would  provide  a  better 
understanding  of  the  plume  rise  and  looping  phenomenon. 

These  data  should  further  enhance  the  predicting  capability 
of  AQAM. 
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